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Abstract 

A model for the Universe is proposed whose constituents are: (a) a dark energy field mod- 
eled by a fermionic field non-minimally coupled with the gravitational field, (b) a matter field 
which consists of pressureless baryonic and dark matter fields and (c) a field which represents 
the radiation and the neutrinos. The coupled system of Dirac's equations and Einstein field 
equations is solved numerically by considering a spatially flat homogeneous and isotropic Uni- 
verse. It is shown that the proposed model can reproduce the expected red-shift behaviors of 
the deceleration parameter, of the density parameters of each constituent and of the luminosity 
distance. Furthermore, for small values of the red-shift the constant which couples the fermionic 
and gravitational fields has a remarkable influence on the density and deceleration parameters. 

1 Introduction 

In the last decade several works appeared in the literature with the aim the explanation of the 
astronomical observations [T] that the Universe is presently expanding with a positive acceleration. 
The great number of the works admits the existence of a component of the Universe - the so-called 
dark energy - which is the responsible for the present acceleration of the Universe, but whose nature 
still remains unknown. The most natural candidate for dark energy was initially a cosmological 
constant (see e.g. [2]), but other models for dark energy were also proposed. The most common 
models, also called quintessence models, make use of scalar fields and barotropic equations of state 
where the ratio between the pressure and energy density assumes negative values (e.g. [3] and the 
references therein) . Other models consider that the dark energy field is described by exotic equations 
of state like the Chaplygin fluid U or the van der Waals fluid [5] . 

Some cosmological models were also investigated in the literature by considering fermions as 
sources of the gravitational field [5]. Recently, the authors [7] proposed a cosmological model in a 
dissipative Universe and showed that the fermions could be the responsible for accelerated regimes 
in the early Universe by simulating an inflaton field and in the present epoch by simulating a dark 
energy field. Another paper on this subject was published very recently [5]. 

In most of the works which investigate the role of the scalar field as source of the gravitational 
field it is common to consider it to be minimally coupled to the gravitational field. However, there 
exist some works in the literature [3] where it was investigated the effects of non-minimally coupled 
scalar fields during the evolution of the Universe. 

The aim of the present work is to describe a spatially flat, homogeneous and isotropic Universe 
characterized by the Robertson- Walker metric and whose constituents are a dark energy field, a 
matter field and a field which represents the radiation and the neutrinos - from now on called 
radiation field. The dark energy field is modeled by a fermionic field which is non-minimally coupled 
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with the gravitational field whereas the matter field represents the pressureless baryonic and dark 
matter fields. The radiation and matter fields are supposed to behave as non-interacting fields. The 
basic field equations for the evolution of the Universe follow from Dirac's equations and Einstein's 
field equations which are solved numerically for given initial conditions. The red-shift behaviors 
of the deceleration parameter, of the density parameters of each constituent and of the luminosity 
distance, which result from the numerical solutions, are compared with the available data set in 
order to verify the viability of the proposed model. 

The manuscript is structured as follows: in Sec. II the Dirac equations and Einstein field 
equations are derived from the action of a non-minimally coupled fcrmionic field. The field equations 
for a flat; homogeneous and isotropic Universe are given in Sec. Ill, which are rewritten in terms of 
the red-shift in Sec. IV. In Sec. V the cosmological solutions are obtained by solving numerically 
the system of field equations for given initial conditions. We close the work with Sec. VI where we 
address to some remarks and sum up the results obtained. The metric signature used is (+, — , — , — ) 
and units have been chosen so that SttG ~ c = h = I. 



2 Non-minimally coupled fermionic field 

We are interested in investigate a Universe modeled by a mixture whose constituents are the fields 
of dark energy, matter and radiation. As was pointed out in the previous section, the dark energy is 
described by a fcrmionic field which is non-minimally coupled with the gravitational field, whereas 
the fields of matter and radiation are considered as non-interacting. The action for this model reads 



i(l-^VV^)i? + £/ + £,„ + £, 



(1) 



where Cm and Cr denote the Lagrangian densities of matter and radiation, respectively. Further- 
more, and ijj are the spinor field and its adjoint, respectively, and ^ is the coupling constant 
between the spinor field and the curvature scalar R. The spinor fields are treated here as classi- 
cally commuting fields [7l [10] . Finally, Cf is the Lagrangian density of the fermionic field which for 
massless fermions reads 

Cf = [v^F^i^^v - (^/.V^)r^V'] - V. (2) 

In cq. ([2]) - which is only a function of and ip - describes the potential density of self-interaction 
between fermions. Moreover, F^ ~ e^j'^ are the generalized Dirac-Pauli matrices, e° denotes the 
tetrad or "vierbein" while the covariant derivatives are given by 

D^ip ^ d^ip ~ n^j,^, Dfj,ijj ^ dfj,ij + ipilfj,. (3) 

Above, the spin connection is defined by 

^^.--J5p.[^^^-e^9,e^]F-F^ (4) 

with denoting the Christoffel symbols. 

From Euler-Lagrange equations applied to the total Lagrangian density in eq. ([T]) it follows the 
Dirac equations for the spinor field and its adjoint coupled to the gravitational field, namely 

iTf^D^ij _ ^ _ = 0, (5) 
dtp ^ 

zi^^V^F^-F^-t-|i?V^=0. (6) 
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The variation of the action ^ with respect to the tetrad leads to Einstein's field equations 

R,.-lg,.R^-j^, (7) 

where T^^ is the total energy-momentum tensor of the sources which is a sum of the energy- 
momentum tensor of the fields of fermions T^"', matter T^'^ and radiation T/"^, i.e., T'"' = T^'^ + 

m ' r • 

The symmetric form of the energy-momentum tensor of the fermions which follows from the 
variation of the action with respect to the tetrad reads 



f ~2 



-CCr^" ~'r\9n- (8) 



Above, the parentheses around the indices denote the symmetric part of a tensor and it was intro- 
duced the tensor T^^, defined by 

3 Field equations for Robertson- Walker metric 

From now on we shall investigate cosmological solutions which can be obtained from the above 
proposed model when applied to a homogeneous, isotropic and spatially flat Universe described by 
the Robertson- Walker metric ds^ ~ dt^ ~ a{t)^5ijdx^dx^ , where a{t) denotes the cosmic scale factor. 

For the Robertson- Walker metric the components of the tetrad, spin connection and Dirac-Pauli 
matrices reduce to 

4 = 6i;, < = ^'5r, Oo = 0, 0, = ia(t)7V, (10) 

r° = 7", = 4t7', = -tV^r^r^r^r^ = 7^ (11) 

Moreover, in this case the spinor field is an exclusive function of time and the Dirac equations ([5|) 
and ([6]) become 

^ + ^iJ^ + ,/^ + ^|i?/^ = 0, (12) 
^ d%! ^ 

?+^i?^-*^7°-4^V^7° = 0, (13) 

where H = a{t)/a{t) denotes the Hubble parameter and the dot represents a differentiation with 
respect to time. 

For a homogeneous and isotropic Universe the total energy-momentum tensor of the sources of 
the gravitational field is written as [Tj^) = diag (p, — p, — p, — p), where p and p denote the total energy 
density and pressure of the sources, respectively. In terms of the energy densities and pressures of 
the constituents, the total energy and pressure are written as 

P = p/+Pm+Pr, P^Pf +Pm+Pr, (14) 

i.e., as a sum of the energy densities and pressures of the fermionic, matter and radiation fields, 
respectively. 
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From Einstein's field equations (O it follows the modified forms of Fricdmann and acceleration 
equations 

g2_ P _ a_ P + '^P 

3(1 -CVV')' a 6(1 -^VV-)' 

respectively. 

By differentiating the Fricdmann equation (jlSP i with respect to time and taking into account: 
(i) the acceleration equation p^ 9. (ii) the Dirac equations p2|) and p3)) and (iii) that the self- 
interaction potential V between the fermions does depend only on the bilinear 4" = tpip, it follows 
the evolution equation for the total energy density of the sources of the gravitational field, namely 

p + 3H{p + p) = -j^^^^. (16) 

The above equation follows also from the Bianchi identities applied to the Einstein field equation 
([7]). Note that due to the coupling between the fcrmionic and the gravitational fields, the right-hand 
side of p6p is not zero and represents the energy transfer between them. As it is well known, only 
two among the three equations given in (jl5p and (|16p arc linearly independent. 

By assuming that the energy-momentum tensors of the constituents have the same form as the 
one for the total energy-momentum tensor, it follows from ([8]) and ([9]) that the energy density and 
pressure of the fcrmionic field read 

Pf = V + 3^H^, (17) 

Note that the existence of a coupling between the fcrmionic and the gravitational fields leads to 
a dependence of the energy density and pressure of the fcrmionic field with respect to the energy 
density and pressure of all constituents through the Hubble parameter H and the scalar curvature 
R. 

After the matter and radiation decoupling both fields behave as non-interacting fields so that 
the evolution equations for their energy densities read 

p„, + 3Hp^ = 0, pr+ iHpr = 0, (19) 

by taking into account that matter is considered as a pressureless fluid, i.e., = whereas the 
radiation obeys the barotropic equation of state Pr = Pr/i- 

The evolution equation for the energy density of the fcrmionic field follows from p^ . pB]) and 
(fig)) , yielding 

pf + m{pf+pf) = -j^^^y (20) 

where the right-hand side of the above equation represents the energy transfer of the fcrmionic field 
to the gravitational field due to their coupling. 

4 Field Equations in terms of the red-shift 

In order to get cosmological solutions for the proposed model, the red-shift will be used as a variable 
instead of time thanks to the following relationships 

z = --l, ± = -H{l + z)^. (21) 
a at dz 
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As a first step we note that that eqs. can be easily integrated leading to the well-known 

dependence of the energy densities of the matter and radiation fields with the red-shift 



Pr{z) = Prim + Z)\ 



(22) 



where Pm(0) and /Or(0) arc the values of the energy densities at present time 2 = 0. 

Afterwards, the energy density and the pressure of the fcrmionic field arc expressed in terms of 
the red-shift, yielding 

Pf^V^'AiH\l + z)^', (23) 
n ~1 



Pf 



2+.e^ 

2(1 -^*) 



dV_ 



(1 + z)H^^' + (1 + zfHH'^'] 



^[(1 + 

e*(/0m+p/) 



(24) 



2(1-^*) 

where the prime denotes the differentiation with respect to the red-shift z and again by assuming 
that the self- interaction potential V between the fermions does depend only on the bilinear ^' = 
For a non-coupled fcrmionic field its energy density and pressure reduce to 



dV 



* - 1/. 



(25) 



Finally, a system of two coupled differential equations for the bilinear ^E* and for the Hubble pa- 
rameter H is obtained from the Dirac equations (|12p , (I13p and from the Fricdmann and acceleration 
equations (|15p . that reads 



(1 + z)H'^" - [2H - (1 + z)H']^' - 3H'^ = 0, 
2iJ(l + z)(l-^*)F' = p + p. 



(26) 
(27) 



For the determination of the cosmological solutions from the above equations we have also to 
specify how the self- interaction potential between the fermions V depend on the bilinear ^. Here 
we suppose that it is given by the power law V = k/[^'(z)]", where k and a are two free parameters. 



5 Cosmological solutions 

For the search of numerical solutions of the system of differential equations (pS]) and ((TT]) together 
with ([22]) through ((24l) we have to specify initial conditions for ^(z), ^''(z) and H{z). Here we shall 
suppose that in the radiation dominated era at zq = 10^ (say) the value of the bilinear is very small 
^(zq) = 10^^ with a smaller value for its derivative with respect to the red-shift ^''(zo) = 10^^'^ 
so that the fcrmionic field rolls very smoothly during the evolution of the Universe. Moreover, 
we shall adopt the following values for the density parameters VLi{z) = pi{z)/p(z) at present time: 
nf{0) = 0.72, nm{0) = 0.279916 and 17^(0) = 8.4 x IQ-^ (see e.g. [TT]). Hence the initial value for 
the dimensionless Hubble parameter can be obtained from Fricdmann equation (jl5|l i together with 
(I22|) and (123|), yielding 



„( / »n^(0)(l + z)3 + r!.(0)(l + z)4 +7c/[vl/(zo)]" 

"^^'"^"V 3(1 -e*(^o)) + 3^(1 + z„)*'(zo) ' ^'^^ 

where k is also a dimensionless quantity. 

To find numerical solutions of the system of differential equations ((26|) and ((27|) , still remains to 
fix values for the free parameters: (a) k and a of the self-interaction potential between the fermions 
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Figure 1 : Density parameters as functions of the red-shift for the ACDM model and for the fcrmionic 
model with different values of Large frame < z < 2, small frame 1000 < z < 4000. 



and (b) the coupling constant between the fermionic field and the gravitational field The following 
values K = 0.6 and a = 1/10 were found by imposing that the present value of the density parameter 
of the fermionic field matches the expected value 5^/(0) = 0.72 when there is no coupling, i.e., when 
5 = 0. The search for values of these parameters by changing the initial conditions is exhaustible, 
specially due to the fact that the system of differential equations ([26|) and ([27|) is very unstable. 
Here we shall adopt the above quoted values. 

In figure 1 the density parameters are plotted as functions of the red-shift in the range between 
< z < 2 (large frame) and in the range between 1000 < z < 4000 (small frame). We infer from 
the large frame of this figure that the coupling parameter ^ has influence on the decay of the dark 
energy field and the corresponding increase of the matter field. Moreover, the density parameter 
of the field of radiation in this range is very small. When ^ = there is no sensible difference 
between the ACDM model and the model where the fermions represent the dark energy field and 
there is no coupling. We note that by increasing the coupling parameter ^ the decrease of the 
fermionic field and the corresponding increase of the matter field are less accentuated. This fact can 
be understood by noting that according to eq. (|20p a part of the energy density of the fermionic 
field is transferred to the gravitational field due to its coupling. Furthermore, one can observe that 
the system of differential equations is very sensible to small changes of the coupling parameter. The 
density parameter of the fermionic field becomes negligible when z « 5 for the case when there is 
no coupling (^ = 0) and when z « 200 for the other two cases ^ = 0.0180 and ^ = 0.0195, with 
a more accentuated decay of the former with respect to the latter. The density parameters of the 
radiation and matter fields are represented in the small frame in figure 1, which shows that the 
equality between the matter and radiation occur at z « 3300. It is noteworthy that here we have 
taken into account the degrees of freedom of the relativistic neutrinos. 

The deceleration parameter q = 1/2 -I- 3p/2p is plotted in figure 2 as function of the red-shift. 
One can infer that the coupling parameter ^ has also a very remarkable influence on the deceleration 
parameter q. It is important to note that the ACDM model furnishes a different solution with 
respect to the uncoupled fermionic model (^ = 0). Indeed, although the values of q coincide at large 
values of the red-shift (z > 1.5) their present values at z = and red-shifts zj that correspond to 
the transition from a decelerated to an accelerated regime read: (a) ACDM model q{0) ~ —0.58 
and Zt = 0.73 and (b) imcoupled model q{Q) ~ —0.68 and zt = 0.80. The corresponding values for 
5^0 are: (c) q{0) = -0.58 and zt = 0.64 for ^ = 0.0180 and (d) g(0) = -0.49 and z* = 0.50 for 
5 = 0.0195. Apart the last case all others fit more or less the experimental values for the deceleration 
parameter, namely, q{0) = -0.74 ±0.18 (see [l2]) and z* = 0.46 ±0.13 (see [T3]). 

In figure 3 it is plotted the difference between the apparent m and the absolute magnitude M of 
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Figure 2: Deceleration parameters as functions of the red-shift for the ACDM model and for the 
fcrmionic model with different values of ^. 
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Figure 3: /xq as function of the red-shift z for the ACDM model and for the fermionic model. Large 
frame < z < 1, small frame < 2: < 0.5. 



a source, denoted by /io and defined by 

Mo = m-il/ = 51og|(l + z)^ _^|+ 25, (29) 

where the quantity between braces represents the luminosity distance in Mpc. The circles in this 
figure correspond to the experimental data for super-novae of type la taken from the recent work by 
Riess et al. [13] . It is important to call attention to the fact that the values of the coupling constant 
^ chosen above do not lead to sensible differences in the values of ^0 so that only one curve for the 
fermionic model is represented in the figure. One can infer from the small frame of figure 3 that 
there is a very good agreement of the curves of the ACDM model and of the fermionic model with 
the experimental data for red-shifts in the range < z < 0.5. From the large frame of this figure one 
observes that there exist some differences between the curves of the two models at higher red-shifts 
0.5 < z < 1. 

6 Final remarks and conclusions 

In this work we have considered a classical fermionic field non-minimally coupled with the gravita- 
tional field. Classical spinors were discussed by Armcndariz-Picon and Greene in Ref. [6] and for 
more details one is referred to this work. Here we point out that: (i) the expectation value of a 
spinorial field in a physical state is a complex number and not a Grassmannian number and (ii) the 
spinorial field can be treated classically if its state is close to the vacuum. 
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Recently, exact cosmological solutions for a Universe described by a Robertson- Walker metric 
and modeled by a neutral scalar field and neutrinos were obtained in the work [15| . The model of 
this work is different from the model presented here since in the former the neutrinos are considered 
as a Majorana fermion field while the neutral scalar field represents the acceleron of the mass- varying 
neutrinos. 

One important point which deserves attention is the one concerning the equivalence of the above 
description with the one in which the Einstein's field equations is written as 



1 



where = + T,^,'' + T/"', with 



(30) 



(31) 



In this case we have Ti^'^.^ = thanks to the Bianchi identities. If we write the total energy- 
momentum tensor as (T^) = diag {p, — p, —p, —p) - so that p and p represent the total energy density 
and the total pressure, respectively - it follows 



P 
3' 



a 
a 



3p 



6 



p + 3H{p + p) ^ 0, 



(32) 



which are the usual form of the Friedmann, acceleration and energy conservation equations, respec- 
tively. 

Equations ([32]) are not in contradiction with the corresponding eqs. (flS]) and (flBl) . since one can 
prove that the energy density and pressure of the latter case are related with those of former case 

by 

Hence, it is straightforward to obtain eqs. (fT5|) and (fT6|l from eqs. (|32|) and ([33|) by the use of 
Dirac's equations and (|13p . 

Another remark is that the non-minimally coupling is important only at small red-shifts and its 
effect becomes negligible at high red-shifts where the matter-radiation decouples and the nucleosyn- 
thesis begins. Moreover, the method adopted here was to solve the system of differential equations 
([26|) and ([27]) by specifying small values for '^{zq) and ^''(^o) at a high red-shift zq so that the 
fermionic field rolls very smoothly during the evolution of the Universe. Another method follows by 
combining equations (fT^ and p3p so that the resulting equation leads to 



*(z) _ (l + z)3 



^(zo) (l + zo)^' 

From the above equation one can obtain that the energy density of the fermionic field reads 



Pf 



K 2.£,[p,n{z)+pr{z)]^{z) 

[vi'(z)]" [i-e*(^)] 



l~^^{z) 



l + 2e*(z) 



(34) 



(35) 



Now by specifying an initial condition to the bilinear '^{z) one can made a similar analysis on the 
behavior of the density parameters, of the deceleration parameter and of the luminosity distance as 
was done in Sec. V. 
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As a last remark we call attention to the fact that we have considered a very simple type 
self-interaction potential which depends only on the bilinear ipip. With this choice, the fermionic 
field seems to behave as a bosonic field, but a comparison of the field equations and solutions of 
the present work with those obtained by Binder and Kremer in Ref. [S] for a bosonic field non- 
minimally coupled with the gravitational field shows that the solutions are quite different. A more 
general self-interaction potential should also depend on the pseudo-scalar invariant (iipj^ip)'^ as the 
one investigated by the authors in Ref. 7: . This work is in progress and shall appear in a forthcoming 
paper. 

To sum up: (a) the proposed model is very sensible concerning the choice of the initial conditions 
and the choice of the free parameters ^, k and a; (b) the density parameter of the fermionic and 
matter fields decrease more slowly with respect to the red-shift when the fcrmions are coupled 
7^ 0) due to a more efficient energy transfer between the fermionic and the gravitational fields; 
(c) by increasing the coupling constant ^ the transition from a decelerated to an accelerated regime 
occurs at earlier red-shifts and (d) there is no significant difference in parameter ^ which is related 
to the luminosity distance - for coupled and uncoupled fermionic fields. 
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